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Preface 


This  work  began  as  an  investigation  into  the  various 
techniques  of  propagation  of  EM  fields  in  optical  devices. 

The  two  simple  techniques,  geometric  optics  and  scalar 
diffraction  theory,  are  both  based  on  assumptions  that 
violate  Maxwell’s  equations,  except  for  special  cases.  The 
current  vector  diffraction  theories  were  rigorous,  but 
difficult  to  apply  to  various  problems.  Major  Glenn  Doughty, 
who  acted  as  my  advisor,  had  the  initial  insight  of  using 
complex  field  vectors  in  the  rigorous  diffraction  problem, 
which  led  directly  to  the  development  detailed  in  the  first 
two  sections  of  Chapter  III.  This  insight  was  so  fruitful 
that  we  decided  to  redirect  the  effort  toward  developing 
the  complex  field  vector  approach  to  diffraction  theory. 

The  significance  of  the  approach  lies  in  (1)  specifying  the 
resultant  field  in  terms  of  separated  component  equations,  and 
(2)  deriving  the  Rayleigh-Sommerfeld  equation.  The  second 
point  helps  explain  where  scalar  diffraction  theory  can  be 
an  accurate  description  of  a  vector  field. 

I  am  most  indebted  to  Major  Doughty  for  his  insight, 
guidance  and  patience  in  our  many  hours  of  discussion.  My 
thanks  is  also  extended  to  Ms.  Sharon  Gabriel  for  her  typing 
assistance.  I  cannot  add  sufficient  appreciation  for  my 
loving  wife,  Jan,  who  supported  me  equally  in  times  of  victory 
and  discouragement.  Most  of  all,  I  thank  my  Lord  Jesus  Christ 
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for  His  strength  and  love.  As  the  Great  Physicist,  He 
has  shown  me  ''great  and  mighty  things"  which  I  did  not 
know.  (Jer  33:3)  I  offer  this  work  to  His  praise  and 
glory. 


Mark  E.  Rogers 
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Abstract 

A  rigorous  solution  to  the  diffraction  problem  is  obtained 
by  using  two  complex  field  vectors:  Q  *>  yH  +  i/yT  E  and 
P  ■  y  M  -  i/yT  E  .  The  derivation  begins  with  Maxwell’s 
equations  that  are  uncoupled  in  terms  of  Q  and  P  .  Then, 
following  the  work  of  Stratton  and  Chu  and  that  of  Doughty, 
the  field  equations  are  integrated  directly  to  yield  a  pair 
of  uncoupled  vector  integral  equations  involving  the  tangential 
components  of  Q  and  P  on  some  open  surface  S'  .  When  S' 
is  planar,  the  equations  are  expanded  in  a  more  useable  set  of 
six  component  integral  equations.  The  main  assumptions  in  the 
derivation  of  these  latter  equations  are  that  the  initial  E 
and  H  must  satisfy  Maxwell's  equation  on  S'  and  that  the 
resultant  field  must  be  calculated  for  k  >>  1/R  .  For  the 

case  where  (1)  the  initial  field  is  traveling  normal  to  S'  , 
and  (2)  the  resultant  field  is  calculated  only  in  the  region 
near  the  optical  axis,  the  expressions  for  the  tangential 
components  of  the  resultant  E  are  identical  in  form  with  the 
Rayleigh-Sommerfeld  equation  of  scalar  diffraction  theory. 

Three  examples  are  developed  to  show  the  method  of  applying 
complex  field  vectors  to  diffraction  problems  and  to  show  the 
agreement  on  axis  with  the  results  calculated  from  the  Rayleigh- 
Sommerfeld  equation.  The  examples  also  provide  insight  into 
the  diffraction  process  by  discussing  the  approximations 
made  in  obtaining  the  resultant  field  near  the  axis.  Thus, 
the  complex  field  vector  approach  provides  a  rigorous,  yet 
simple  and  straightforward,  method  of  solving  the  diffraction 
problem. 
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THE  USE  OF  COMPLEX  FIELD  VECTORS 
IN  DIFFRACTION  THEORY 

I .  Introduction 

The  goal  of  the  general  diffraction  problem  is  to 
calculate  a  resultant  electromagnetic  field  on  some  surface 
S  from  an  initial  electromagnetic  field  on  some  surface 
S'  .  Building  on  the  work  of  Huygens,  Fresnel,  Kirchoff 
and  others,  many  authors  (e.g.,  Refs  1-5)  have  taken  various 
approaches  to  solving  this  problem.  However,  the  two  basic 
approaches  are  both  based  in  principle  on  Maxwell's 
equations.  The  first  approach  deals  rigorously  with  the 
vector  nature  of  the  EM  field,  and  is  known  as  vector 
diffraction  theory.  The  second  approach  assumes  that  each 
field  component  may  be  considered  independently.  This 
scalar  diffraction  theory  is  suggested  by  the  separate  wave 
equations  for  the  electric  and  magnetic  field  intensities, 

F  and  H  ,  that  are  obtainable  from  Maxwell's  equations. 
However,  the  "real  problem  is  not  the  integration  of  a  wave 
equation,  either  scalar  or  vector,  but  of  a  simultaneous 
system  of  first-order  vector  equations  relating  the  vectors 
E  and  H  ."  (Ref  6:463)  Unfortunately,  the  inherent 
coupling  of  E  and  H  makes  the  rigorous  solution  of  the 
general  diffraction  problem  difficult  to  obtain. 
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The  vector  diffraction  theories  have  some  common 
advantages  and  disadvantages.  Each  is  a  rigorous  solution 
to  Maxwell's  equations  for  a  given  set  of  assumptions  and 
boundary  conditions.  However,  the  resulting  expressions 
are  often  difficult  to  apply  to  general  problems  and  may 
obscure  the  underlying  physical  phenomena  involved  in 
diffraction. 

Scalar  diffraction  theory  is  a  more  phenomenological 
approach,  but  it  is  more  easily  applied  and  quite  accurate 
in  certain  regions.  The  approach,  which  is  a  rigorous 
description  of  the  propagation  of  sound  waves,  is  suggested 
by  the  separation  of  the  vector  wave  equations  for  E  and 
H  into  six  scalar  wave  equations  when  rectangular  coordinates 
are  used.  However,  the  assumption  that  each  component  is 
independent  ignores  the  inherent  coupling  of  E  and  H  , 
and  the  results  generally  do  not  satisfy  Maxwell's  equations. 

From  the  above  discussion,  we  see  that  no  approach  is 
both  rigorous  and  easily  applied.  The  goal  of  this  work  is 
to  develop  a  theory  of  diffraction  that  incorporates  the 
advantages  of  the  vector  and  scalar  approaches  while  avoiding 
some  of  the  shortcomings.  Specifically,  the  primary  objective 
of  this  work  is  to  develop  an  approach  to  diffraction  that 
uses  complex  field  vectors  to  recast  Maxwell's  equations 
into  a  set  of  separated  integral  equations  that  describe  the 
resultant  EM  field.  Two  secondary  objectives  are:  (1)  to  show 
that  scalar  diffraction  theory  is  a  special  case  of  this 
approach,  and  (2)  to  show  that  this  approach  is  readily 
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applicable  to  various  diffraction  problems.  A  minor 
objective  in  this  application  is  to  gain  some  familiarity 
with  Fast  Fourier  Transform  techniques  as  applied  to  the 
computation  of  solutions  to  diffraction  problems.  The 
objective  of  this  report  is  to  describe  the  development  and 
application  of  this  approach  to  the  general  diffraction  problem. 

The  use  of  complex  field  vectors,  Q  =  pH  +  i/y"e  E  and 
P  *  yH  -  i/yT  E  ,  has  two  significant  features.  The  first 
feature  is  the  development  of  a  set  of  uncoupled  vector 
integral  equations  that  retain  the  vector  nature  of  the  EM 
field.  The  component  equations  that  result  from  this 
decoupling  are  easily  applied,  interpreted  and  understood, 
and  yet  they  are  still  rigorously  based  on  Maxwell's  equations. 
The  second  contribution  is  the  substantiation  of  scalar 
diffraction  theory  under  conditions  common  to  many  diffraction 
problems. 

This  presentation  begins  with  a  review  of  past  work  and 
a  discussion  of  some  theoretical  preliminaries,  both  of 
which  are  contained  in  Section  II.  This  review  is  intended 
to  set  the  stage  for  later  developments,  and  is  not  an 
exhaustive  history  of  diffraction  theory.  The  theoretical 
preliminaries  provide  the  background  for  the  theoretical 
development  of  the  complex  field  vector  approach  in  Section 
III.  Three  example  fields  are  examined  in  Section  IV  to 


illustrate  the  application  of  this  approach.  Section  V 
contains  the  conclusions  and  recommendations.  Various 
appendices  contain  detailed  mathematical  developments  to 
support  this  work.  All  equations  will  be  in  MKS  units 
unless  otherwise  specified. 
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II.  Theoretical  Preliminaries  and  Past  Work 


This  section  develops  the  geometry  for  the  general 
diffraction  problem,  discusses  the  type  of  medium  assumed  in 
this  report,  and  derives  some  preliminary  relations  based 
on  Maxwell's  equations.  The  section  also  includes  a 
discussion  of  two  areas  of  past  work  that  are  relevant  to 
this  report.  The  first  is  a  direct  integration  of  Maxwell's 
equations  as  developed  by  Stratton  and  Chu  (Ref  4)  and 
extended  by  Doughty  (Ref  7);  the  second  is  the  Rayleigh- 
Sommerfeld  formulation  of  scalar  diffraction  theory. 

The  reader  interested  in  a  more  in-depth  treatment  of 
vector  diffraction  theory  is  referred  to  the  paper  by 
Kottler  (Ref  2),  the  article  by  Bouwkamp  (Ref  S),  the 
treatment  by  Born  and  Wolf  (Ref  8:556),  or  the  book  by 
Baker  and  Copson  (Ref  9:102).  Scalar  diffraction  theory  is 
discussed  in  most  texts  on  physical  optics,  such  as  the  text 
by  Born  and  Wolf  (Ref  8:320).  The  Rayleigh-Sommerfeld 
formulation  described  in  this  section  is  based  on  the 
discussion  in  Goodman  (Ref  10:30). 

Theoretical  Preliminaries 

Figure  1  shows  the  geometry  for  the  general  diffraction 
problem  used  in  this  report.  The  initial  EM  field  is  assumed 
to  be  non-zero  only  on  the  open  surface  S"  ; 
values  on  S'  will  be  denoted  using  primes.  The  resultant 
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Figure  1.  Geometry  for  the  General  Diffraction  Problem 


EM  field  is  described  on  some  open  surface  S  ;  values  on 
S  will  be  unprimed.  The  distance  from  a  source  point  on 
S'  to  an  observation  point  on  S  is  given  by  R  ■  | r - F'  | 
and  the  angle  e  is  defined  by  the  angle  between  r ~ r ' 

A  A 

and  n  ,  where  n  is  a  unit  normal  to  S'  pointing 

toward  the  surface  S  .  The  origin  is  chosen  to  lie  on 

the  S'  such  that  if  S'  and  S  were  plane  surfaces,  the 
points  on  S'  would  be  described  by  (x',y',0)  and  those 
on  S  by  (x,y,z)  .  For  many  applications,  S'  and  S 

will  be  planes,  but  in  the  general  equations  they  need  not 

be  so  restricted. 

Throughout  this  report,  the  propagation  medium  will  be 
assumed  to  be  homogeneous,  isotropic,  linear  and  free  of 
charge,  unless  otherwise  stated.  Ohm’s  Law  is  assumed  to 
hold  in  the  medium,  unless  otherwise  stated.  Under  these 
assumptions.  Maxwell’s  equations  can  be  written  as  follows: 


7  x  E  (r,t)  -  -y  £  [H(r,t)]  (1) 

V  x  H  (r,t)  -  e  £  [E(r ,t)]  ♦  oE  (2) 

V  •  E  (r,t)  •  0  (5) 

yV  .  IT  (r,t)  •  0  (4) 
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Both  I"  and  H  can  be  expanded  here  in  normal  modes,  as 
shown  in  Eq  (5)  for  an  arbitrary  vector  A  . 

Oft 

A(r  ,t)  -  — - —  /  A(r,u)  elu>t  du  (5) 

(2n)  3  -«• 

Using  these  expansions  in  Eqs  (1)  and  (2),  one  finds  for  an 
arbitrary  frequency  u  , 


v  x  E(r) 

+  ioiy 

H(r)  -  0 

(6) 

v  x  H(r) 

-  io>£ 

E(r)  -  0 

(7) 

where  the  frequency  dependence  has  been  suppressed  and  the 
complex  dielectric  constant  is  defined  as  e  *  e  -  i  —  . 

ii> 

Equations  (6)  and  (7)  provide  the  starting  point  for  the 
development  in  Section  III.  Before  this  development  begins, 
however,  it  will  be  useful  to  review  some  past  work  that  is 
relevant  to  this  report. 

Past  Work 

A  Vector  Diffraction  Theory.  Stratton  and  Chu 
developed  an  analytic  solution  to  the  diffraction  problem 
by  direct  integration  of  Maxwell’s  equations  using  a 
vector  analogue  to  Green's  theorem  (Ref  4  or  Ref  6:464). 

! 

For  a  nonconducting  medium,  they  obtained  the  resultant 
E  and  H  at  any  point  on  the  interior  of  a  volume  V 
bounded  by  a  closed  surface  S£  ,  where 
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E  and  H  are  expressed  as  the  sum  of  a  volume  integral  and 
an  integral  over  Sc  .  The  initial  field  is  specified  on 
Sc  .  By  assuming  the  volume  to  be  free  of  charge  and  current 
and  the  initial  field  to  be  non-zero  only  on  some  finite 
portion  of  SQ  ,  say  S'  ,  the  integral  equations  for  E 
and  R  can  be  rewritten  as  the  sum  of  an  integral  over  S' 
and  an  integral  on  the  contour  c'  enclosing  S'  .  These 
equations,  Eqs  (8)  and  (9),  are  equivalent  to  those  derived 
by  Kottler  (Ref  2). 

E(r)  -  -r-  iv*.  H'.  d£ 

*  luc  4  IT  1*1 


1 

*7 


/  [iuu  (nxH')*. 
S'  1 


♦ 


A 

(nxl')  x  + 


(n-E^V^ldS' 

(8) 


H(r)  - 


1 

luu 


*  d? 


♦ 


[iue (r  x 


(n  x  H')  x 


(n*ff')v*1]dS' 


(9) 


where  »  exp  [ikR]/R  .  In  a  parallel  development, 
Doughty  (Ref  7:21)  developed  an  analogous  set  of  equations 
for  the  case  where  Ohm’s  Law  holds  in  the  medium  and  without 
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assuming  the  initial  field  is  zero  everywhere  except  on  S'. 
For  the  diffraction  problem,  Doughty's  approach  is  equivalent 
to  that  of  Stratton  and  Chu.  Taking  the  curl  of  Doughty's 
equations,  one  obtains 

E(r)  »  f  (nxH")  •vxv'xC  +  vxG •  (nxE  )]dS'  (10) 

S 


H(r) 


[^(nxE')*vxv'xS‘  +  vxC- (nxH')]dS' 


(11) 


where  is  a  complex  dielectric  constant  that  accounts 

for  the  absorption  or  gain  in  a  medium,  and  v'  refers  to 
the  gradient  with  respect  to  the  primed  coordinates 
(x',y',z')  ,and  G  is  the  dyadic  Green’s  function. 


T  • 


(a  a  + 
xx 


Vy  *  V»’ 


(12) 


where 


*  *  exp  [-  ikR]/4*R  (13) 

where  R  »  ) r  -  r ' |  .  The  resultant  EM  field  expressed  by 

Eqs  (10)  and  (11)  is  a  solution  to  Maxwell's  equations  but, 
owing  to  the  coupling  of  E  and  H  ,  the  system  of  equations 
is  not  readily  solved. 


i 


10 


Rayleigh-Sommerfeld  Formulation.  A  brief  discussion 
of  scalar  diffraction  theory  highlights  its  phenomenological 
nature,  and  provides  background  for  Sections  III  and  IV.  The 


underlying  assumption  is  that  light  can  be  treated  as  a 
scalar  phenomenon,  described  by  some  scalar  wave  function 
u(r)  that  satisfies  the  scalar  Helmholtz  equation, 

(V2  +  k2)  u(r)  -  0  (14) 


where  the  wave  number  k  -  w/ITe  -  .  u(r)  is 

phenomenologically  interpreted  as  one  of  the  transverse 
components  of  E  or  H  .  By  using  Green's  theorem,  the 
value  of  u(r)  within  a  volume  V  can  be  represented  as  an 
integral  over  a  closed  surface  Sc  that  depends  on  the 
values  of  u  on  S  and  a  scalar  Green's  function  G 
that  satisfies  Eq  (14).  The  closed  surface  can  be  partitioned 
into  a  plane  surface  Sj  and  a  hemispherical  surface  S2 
The  integral  on  S£  approaches  zero  as  the  surface  is  expanded 
toward  infinity,  provided  u(r)  meets  the  Sommerfeld  radiation 
condition.  Then  the  resultant  function  is  represented  by 

U(F)  ■  37/  ucf-)  dSj  os; 

1  * 
where  —  denotes  the  derivative  with  respect  to  n  .  The 

dll 

infinite  plane  is  usually  broken  into  two  regions,  an 

aperture  (such  as  the  open  surface  S'  )  and  the  region 
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exterior  to  the  aperture.  Choosing  G  to  be  $  as  defined 

by  Eq  (13)  leads  to  specifying  both  u  and  —  to  be  zero 

o  ri 

except  on  S'  ,  which  are  inconsistent  boundary  conditions. 
Sommerfeld  showed  that  a  Green's  function  exists  that  is 
zero  on  the  plane  surface  and  thus  the  boundary  conditions 

are  reduced  to  specifying  u(r')  to  be  zero  except  on  S' 

Then  Eq  (15)  becomes 


u  (r) 


i  -ikR 

T  ~ R—  cos 


e  dS' 


(16) 


where  6  was  defined  in  the  preceding  section.  Using 
Eq  (13),  the  above  equation  can  be  rewritten: 


u(r)  -  ik  /  u(r')  (2  cos  0)  *  dS' 
S' 


(17) 


Both  Eqs  (16)  and  (17)  are  forms  of  the  Rayleigh-Sommerfeld 
equation,  which  is  mathematically  consistent,  but  was  obtained 
by  assuming  each  field  component  is  independent.  In  spite 
of  this  deficiency,  the  irradiance  patterns  do  agree  with 
experimental  results  when  e  is  small  (Refs  11,  12). 

Equations  (10),  (11),  and  (17)  are  important  results 
that  will  be  used  in  the  development  and  application  of  the 
complex  field  vector  approach.  This  development  is  the 
subject  of  the  next  section. 
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III.  Theoretical  Development  of  the  Complex 
Field  Vector  Approach  in  Diffraction 

The  purpose  of  this  chapter  is  to  solve  the  diffraction 
problem  with  complex  field  vectors.  In  the  first  section 
below,  these  vectors  are  introduced  and  the  resultant  electro¬ 
magnetic  field  is  expressed  by  a  set  of  uncoupled  vector 
integral  equations  over  open  surfaces.  The  integral  equations 
are  then  specialized  to  open  planar  surfaces;  so  then  the 
resultant  field  is  described  by  six  component  integral 
equations.  A  discussion  of  these  six  equations  provides 
insight  into  the  diffraction  process.  The  chapter  closes  by 
deriving  the  Rayleigh- Sommerf eld  equation  from  these  component 
equations. 

The  Uncoupled  Field  Equations 

As  stated  earlier,  several  different  approaches  could 
be  taken  to  rigorously  solve  the  diffraction  problem.  Since 
most  practical  problems  are  based  on  a  knowledge  of  the  initial 
field,  a  solution  based  directly  on  that  knowledge  seems  most 
useful.  For  this  reason,  the  approach  taken  here  follows 
that  of  ration  and  Chu  (see  Section  II)  and  relies  on  a 
knowledge  of  the  initial  field  intensity  vectors,  E  and  H  . 

Let  the  complex  field  vectors  be  defined  as  in  Eqs  (18) 
and  (19). 
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yH(r)  +  i  /gT  E"(r) 


(18) 


Q(r)  - 

P(r )  -  gH(r)  -  i  /gT  E(r)  (19) 

Others  have  used  combinations  of  E  and  H  to  reduce 
Maxwell’s  equations  to  particularly  compact  form  (Refs  9, 

13,  14).  The  definitions  used  here  follow  Stratton 
(Ref  4:32),  but  the  forms  in  Eqs  (18)  and  (19)  were  first 
suggested  for  application  to  the  diffraction  problem  by 
Doughty  (Ref  15).  Note  that  these  complex  field  vectors 
retain  the  proper  phase  relationship  in  the  field.  In 
P  ,  E  is  advanced  by  90°;  in  Q  ,  E  is  retarded  by 
90°.  Equations  (18)  and  (19)  are  solved  for  E  and 
H  to  obtain  the  equations, 


E  - 

-  %  (q  -  p) 

(20) 

H  - 

17  CQ  ♦  P) 

(21) 

where  v  ■  l//ye  .  Substituting  these  expressions  into 
Eqs  (6)  and  (7),  one  can  solve  the  equations  simultaneously 
to  obtain  uncoupled  equations  in  Q  and  P  : 


Applying  these  equations,  as  well  as  the  definitions  of 
Q  and  P  ,  Eqs  (10)  and  (11)  can  be  manipulated  by  first 
multiplying  Eq  (10)  by  i/yT  and  Eq  (11)  by  y  ,  and  then 
adding  or  subtracting  Eq  (10)  from  Eq  (11)  to  yield  two 
uncoupled  integral  equations  over  open  surfaces  that  depend 
on  the  values  of  Q  and  P  on  the  open  surface,  as  shown 
in  Eqs  (24)  and  (25). 

Q(r)  -  /  {[nxQ'(r')]-VxC  +  ^[nxQ(F')]  •  vxv'xG-}  dS'  (24) 

s 

P(r)  «  / J [nxP '(?')]• vx&  -  £[nxP'(r ')] -vxv'xC}  dS'  (25) 

S 

Recall  that  G  is  the  dyadic  Green’s  function  defined  in 
Eq  (12).  The  above  equations  represent  the  resultant  field 
due  to  the  initial  field  on  S'  ,  since  E  and  H  can 
be  obtained  via  Eqs  (20)  and  (21). 

To  properly  use  Eqs  (24)  and  (25),  the  initial  field 
used  to  form  Q'  and  P'  must  be  a  valid  EM  field  on 
S'  .  Thus,  either  E'  ,  H'  or  the  tangential  components 
of  E'  and  H'  on  the  surface  S'  must  be  known  in  order 


r 


to  obtain  the  resultant  field.  The  same  knowledge  of  the 
initial  field  would  permit  Eqs  (10)  and  (11)  to  be  solved, 
but  Eqs  (24)  and  (25)  are  simpler  to  use  because  they  are 
uncoupled.  But  even  these  equations  are  cumbersome.  In 
the  next  section,  we  will  use  these  equations  to  generate 
a  set  of  six  component  integral  equations  on  S'  that  are 
not  only  quite  easy  to  apply,  but  reveal  a  good  deal  about 
the  diffraction  process. 

Application  to  Plane  Surfaces 

A  A 

Let  S'  be  a  plane  surface,  and  choose  n  «  +a^ 

Consider  first  the  integrand  in  Eq  (24).  By  replacing 

*  A  A  __ 

n  with  +a  ,  the  cross  product  of  n  and  Q  decomposes 

4 

Q  into  transverse  components.  The  continued  expansion  of 
the  integrand  into  components  is  facilitated  by  recalling 
R  *  |r  -  r ' |  which  leads  to  the  relationship  ■  -V'* 
Calculating  the  derivatives  and  taking  the  scalar  product 

_  __  A  A  A 

of  Q(r)  with  a  ,  a  and  a  sequentially,  Eq  (24) 

x  y  z 

can  be  decomposed  into  three  scalar  equations:  (Appendix  A 
contains  details  of  this  decomposition.) 


<V 


I  (-<);•  z 

S'  • 


C  *  7  9  * 


/  <Qx*  X 

s  *  ** 


-  ^yk2**\*x:>  -  q;  Uv]  1 4S' 

-  ♦  V.yx”  dS' 

*  V.y  *  r[Qi‘.xy  •  ds' 


(26) 


(27) 


(28) 
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where  the  notation  on  the  terms  containing  •  represent 
the  derivatives  of  $  with  respect  to  the  subscripted 
variable;  for  example,  t  x  -  fi  and  «>xy  -  . 

Assuming  k  >>  1/R  ,  which  restricts  the  evaluation  of  the 
resultant  EM  field  to  be  at  least  several  wavelengths  from 
surface  S'  ,  the  generalized  first  and  second  derivatives 
are  approximated  by  Eqs  (29)  and  (30),  respectively. 

(See  Appendix  B.) 


x.  -  x: 
ik*  -i.  ,  ■  1 


’xixj 


ai...  - 

R  D  2 


where  x^  is  associated  with  the  values  x  ,  y  ,  and  z 
as  l  ranges  from  1  to  3  ,  respectively.  Using  these 
derivatives  and  the  definitions  below 


sin  a  = 


x  -  x 


;  sin  b  =  £  ^  L  .  ;  cos  e  s  z  ^  z 


one  obtains  three  integral  expressions  for  the  components 
of  Q  ,  as  indicated  in  Eqs  (31)  through  (33) 


Q„  ■  k  /  *[Q'(i  cose-sinasine)  "  Q'(l-sin2a)]  dS'  (31) 

A  g  >  X  j 


0  ■  k  /  ♦ [ Q ' ( i  cose+sinasing)  ♦  Q'(l-sin2g)]  dS' 

/  c  *  y  x 


Q_  ■  k  /  «[-Q'(i  sina+singcose )  -  Q'(i  sing-sinacose)]  dS 
z  S'  x  y 


These  equations  are  written  in  a  particularly  compact 
form  by  introducing  the  complex  components  defined  in 
Eq  (34)  with  an  analogous  definition  for  the  components  P± 
These  components  are  an  alternate  way  of  decomposing 
and  P  into  components.  Although  their  transformation 
properties  are  more  complicated  than  their  rectangular 
counterparts,  these  components  lend  particular  insight 
into  the  diffraction  process,  as  will  be  discussed  in  the 
next  section. 


Q±  *  Qx  »  i  Qy 


C34) 


Using  this  definition  as  well  as  the  definition  that 
6  ■  tan-1  (sina/sine)  and  the  angular  relationship  sin2fi  ■ 
sin2a  +  sin2g  which  is  obtained  from  the  definition  of 
R  ,  Eqs  (31)  through  (33)  become  Eqs  (35)  through  (37). 


ik  /  [Q;(cose+l  -  2iyi*L)  -  Q'  sin2 6  e'2l6]$  dS' 

s 


(35) 


ik  /  [Q;(cose-l  ♦  ^y-)  ♦  jQ;  sin2e  e2l<5]*  dS' 

s 


(36) 


-  |  /  [<^(l+cose)ei<5  -  Q_'(l-cose)  e'ld3*  sine  dS' 

s 


(37) 
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By  a  similar  series  of  steps,  Eqs  (38)  through  (40)  can 
be  obtained  for  P(r)  : 


ik  /  [p^(cose-r  + 


ik  /  [P' (cos 0+1  - 
S' 

T  /  [P+ (1-cos e)e 
S' 


£iy^-)  +  y  P_'  sin2e  e'2i6]$  dS' 

^y^)  -  y  p;  sin20  e2i6]«  dS' 

5  -  P'(l+cos0)  e"*{]«  sine  dS' 


(38) 


(39) 


(40) 


Equations  (35)  through  (40)  are  the  six  component  integral 
equations  that  are  used  to  find  the  resultant  EM  field  if 
S'  is  a  plane.  The  initial  E  and  H  are  assumed  to 
satisfy  Maxwell's  equations,  and  the  component  equations 
were  developed  with  the  approximation  k  >>  1/R  .  One 

should  note  here  that  the  methodology  of  using  the  complex 
field  vectors  in  diffraction  begins  by  taking  E  and  H 
on  plane  S'  and  forming  Q'  and  P'  .  These  vectors 
retain  the  coupling  of  E  and  H  ,  but  allow  the  resultant 
field  to  be  expressed  by  the  six  component  equations  i,iven 
above.  Using  the  definitions  of  Q±  and  P±  ,  Eqs  (20) 
and  (21)  can  be  rewritten  in  a  form  that  facilitates  reforming 
the  resultant  field  in  terms  of  E  and  H  .  Eqs  (41) 
through  (43)  show  the  relationships  for  E  . 
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-  X  (Q+  +  Q-  '  p+  ■  p-) 


T  (Q_  -  Q+  -  P.  ♦  P+) 


T  ‘  V 


The  set  of  scalar  integral  equations  is  a  mathematically 
simple  yet  rigorous  solution  to  the  diffraction  problem. 

Discussion  of  Component  Integral  Equations 

In  this  section  we  analyze  three  terms  that  occur  in 
Eqs  (35)  through  (40),  namely  the  complex  components  Q± 
and  P  ,  the  term  6  ,  and  the  angular  dependence  (which 
is  also  called  the  obliquity  factor). 

The  Complex  Components.  To  investigate  the  nature  of 
Q±  and  P  ,  we  discuss  their  definition  using  polar 
coordinates,  their  expression  in  terms  of  impedances,  and 
their  circularly  polarized  nature. 

First,  the  components  can  be  readily  expressed  in 
cylindrical  coordinates,  as  shown  in  Eq  (44): 


q*  ■  cq  •  i  <y  •’ 


where  the  transformation  matrix  in  Eq  (45)  has  been  used. 
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(45) 


The  expression  in  Eq  (44)  and  the  analogous  expression  for 
P  +  are  useful  in  problems  with  cylindrical  symmetry,  since 
the  component  integral  equations  can  still  be  used  in  their 
present  form. 

Next,  the  complex  scalar  Q+  is  rewritten  in  terms 
of  intrinsic  and  wave  impedances,  using  Eqs  (18)  and  (34) 
as  well  as  the  definition  of  the  intrinsic  and  wave 
impedances  (Ref  16:87): 


+  n 


xy+ 


±ti  + 
xy£ 


(47) 


yx' 


+  ^yx£ 


(48) 


21 


The  +  sign  on  the  subscripts  of  the  wave  impedances 
refers  to  positive  z  propagation,  while  the  -  sign 
refers  to  negative  z  propagation.  The  order  of  x  and 
indicates  which  ratio  of  components  is  implied.  Finally, 
the  combined  symbol  using  the  ambiguous  sign  allows  the 
equations  using  the  wave  impedances  to  be  written  in  a 
compact  manner.  The  resulting  equations  for  Q+  and  Q 
are  shown  in  Eqs  (49)  and  (50).  (See  Appendix  C.)  An 
analogous  procedure  can  be  applied  to  P 


fuli  E  (1±-  )  +  i  /ut  E  (lt>-  ) 

y  n  .  x  n  . 


yx 


xy 


(49) 


Q. 


E  (1+^  )  ♦  i  /ne  E  ) 


-  n 


yx 


xy 


(50) 


Recalling  the  time  factor  is  e 


iut 


E  ■  E+a+  +  E_a_  where 


,  and  writing 


a±  •  'ax 1 1  y 


(51) 


one  can  show  that  E+  is  associated  with  right  circular 
polarization  and  E_  with  left  circular  polarization  for 
+a  propagation,  and  vice  versa  for  -z  propagation. 
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For  uniform  plane  waves,  |n  |  *  |n  |  ■  |n|  .  Then, 

xy  yx 

for  a  uniform  plane  wave  propagating  in  the  +z  direction, 
Eqs  (49)  and  (50)  become 


Q  -  2i  /7e  (E  ♦  i  E  )  -  2/1  i  /y7  E_ 

A  Y 


Q. 


(52) 

(53) 


The  expressions  for  P±  are  P+  *  0  and 

P_  -  -2i  /ITe(Ex  •  i  Ey)  «  -2/2  i  /^TT  E+  (54) 

A  similar  analysis  for  the  -z  direction  reveals  that  Q± 
is  identified  as  left  circular  and  P±  as  right  circular. 

For  predominantly  +z  propagation,  where  k  *  k  a  , 

it 

Q_  and  P+  are  negligible,  while  Q+  and  P_  are 
negligible  for  -z  propagation.  Since  the  quantities  Q± 
and  P±  are  components,  they  do  contain  the  vector  sense 
of  direction  and  polarization.  Any  polarization  state  can 
be  formed  with  P±  and  Q±  by  the  proper  choice  of  E+ 
and  E_ 

We  have  seen  that  the  complex  components  can  be  expressed 
in  either  rectangular  or  cylindrical  coordinates.  The  components 
are  dependent  on  both  the  medium  and  the  ratio  of  the  field 
vectors,  and  they  contain  the  vector  nature  of  the  field  as 
shown  by  the  inherent  circular  polarization  and  directionality 
of  Q±  and  P±  for  predominantly  ±z  propagation. 


The  6  Term.  This  term  represents  a  spatial  phase  shift 
between  the  initial  field  and  the  resultant  field  since  <5 
can  be  expressed  as 

«  *  (y-C  fo  CSS) 

This  term  is  the  only  one  in  the  integrand  that  is  linear 
in  (x-x")  and  Cy-y')  .  It  is  also  worth  noting  that  each 

term  containing  6  is  also  directly  related  to  the  angle 

e  ,  so  as  e  gets  small,  the  sense  of  orientation  is 
weakened.  This  explains  in  part  why  scalar  diffraction  theory 
is  accurate  in  the  paraxial  region  (where  9  is  small). 

Obliquity  Factors.  A  comparison  of  the  various 
obliquity  factors  sets  the  stage  for  the  derivation  of  the 
Rayleigh- Sommerfeld  equation  in  the  next  section.  Consider 
again  an  initial  field  traveling  in  the  predominantly  +z 
direction.  Then,  for  angles  less  than  about  45°,  Eqs  (35) 
and  (39)  dominate  Eqs  (36)  and  (38)  in  terms  of  magnitude, 
since  Q'  >>  Q'  and  ?'  >>  ?'  and  cose  >  sin2e 

Notice  that  the  term  (cos0  ♦  1  -  — y - )  is  common  to  both 

Eqs  (35)  and  (39).  This  term  can  be  compared  against  the 
obliquity  factor  of  the  Fresnel-Kirchhoff  formulation  of 
scalar  diffraction  theory,  1  +  cos9  ,  and  that  of  the 
Rayleigh-Sommerfeld  formulation,  2cos8  .  At  first  glance 
the  new  factor,  identified  as  B  in  Figure  2,  seems  to  be 
a  corrected  Fresnel-Kirchhoff  factor,  identified  as  A  in 
Figure  2.  But  in  the  paraxial  approximation, 
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COS0  ~ 


sin  9  ~  6 


the  Rayleigh-Sommerfeld  factor,  designated  C  in  Figure  2, 


is  identical  with  factor  B 


The  closeness  of  the  curves 


for  B  and  C  suggests  that  the  scalar  integral  equations, 
Eqs  (3S)  and  (38),  may  be  reducible  to  the  Rayleigh-Sommerfeld 
equation,  Eq  (17). 

Derivation  of  Rayleigh-Sommerfeld  Equation 


Consider  again  the  case  of  an  initial  field  propagating 
predominantly  in  the  +z  direction.  As  shown  above, 

Qj  and  P'  are  negligible  in  this  case.  Since  curve  B  and 
curve  C  in  Figure  2  are  so  close,  we  use  the  paraxial  approxi¬ 
mation  of  Eq  (56)  and  write  Ex  by  substitution  of  Eq  (35) 
into  Eq  (41): 

E  -  )£/  {(q;-p;)(2-92)  ♦  [q;  e2i5-p:e‘2i6]|2}*  dS'  (57) 

x  4  S' 


Since  the  term  (2-e2)  can  be  replaced  by  2cose  ,  and  a 
more  rigid  angular  restriction  is  made  such  that  the  terms 

a  2 

containing  -y  are  negligible,  then  the  above  equation 
can  be  written  as  Eq  (58). 


ik  /  E'  (2cos9)«  dS' 
S'  x 


OBLIQUITY  FACTORS 
-D.00  0.75  1.50  2.25  3.00 

*~r  l  I  1  I  _ I 


Figure  2.  Plot  of  Various  Obliquity  Factors 
where  A  *  l+cos0  ,  B  ■  A-isin*0  , 
and  C  *  2cos0  . 
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This  is  the  Rayleigh-Sommerfeld  equation,  Eq  (17),  if  the 
scalar  wave  function  u(r)  is  associated  with  Ex  .  A 
similar  equivalence  is  found  for  E  ,  showing  that  the 
Rayleigh-Sommerfeld  equation  can  be  accurately  applied  to 
the  transverse  components  of  the  initial  field,  provided 
that  (1)  the  initial  field  is  valid  EM  field,  (2)  k  >>  1/R  , 

(3)  k  ;  k  ,  and  (4)  the  angle  e  is  small.  This  last 

it 

restriction,  termed  the  near-axial  approximation,  can  be 
stated  as  e  <<  1  rad  .  Yet,  even  this  approximation  is 
not  overly  restrictive  for  most  applications.  This  derivation 
indicates  why  the  scalar  diffraction  theory  is  an  accurate 
representation  for  the  transverse  components  in  the  near 
axial  region.  The  examples  in  Section  IV  will  illustrate  how 
the  complex  field  vector  approach  is  applied  to  three  different 
polarization  states,  as  well  as  illustrate  that  it  yields 
results  identical  to  those  predicted  by  scalar  diffraction 
theory. 
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IV.  Application  of  the  Complex  Field  Vector  Approach 

The  goal  of  this  section  is  to  illustrate  now  useful 
the  complex  field  vector  approach  is  in  solving  diffraction 
problems.  We  begin  with  some  comments  on  methodology  and 
other  topics  common  to  the  three  examples  used  in  this  section. 
The  three  examples  are  then  presented,  beginning  with  a 
linearly  polarized  field  and  followed  by  fields  of  azimuthal 
and  circular  polarizations.  The  section  concludes  with  some 
observations  on  the  application  of  the  complex  field  vector 
approach  to  diffraction  problems. 

Preliminary  Comments 

As  noted  in  Section  III,  the  initial  E  and  H  are 

specified  on  surface  S'  .  In  each  example  in  this  section, 

the  initial  fields  are  chosen  so  that  the  fields  are  zero 

at  the  tdge  of  S'  .  This  boundary  condition  eliminates 

the  need  to  consider  edge  currents.  The  propagation  vector 

of  the  initial  field  is  assumed  to  be  k  a  ka_  ,  so  that 

z  z 

the  field  is  initially  traveling  in  a  predominantly  +z 
direction.  In  each  example,  a  characteristic  scalar,  F  , 
is  obtained  that  is  proportional  to  the  term  (k  ±  k2) 

If  the  area  S'  has  a  minimum  dimension  that  is  at  least 
several  wavelengths  in  size,  then  F_  can  be  neglected  in 
comparison  to  F+  ,  simplifying  the  analysis.  The  first 
example  will  show  this  comparison  in  detail.  Further,  each 
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example  will  be  solved  by  using  the  far  field  approximation 
which  is  valid  when  the  area  of  S'  ,  denoted  as  A2  , 
is  related  to  the  wavelength  and  distance  to  the  observation 
plane,  z  ,  by  Eq  (59). 


The  area  A2  is  interpreted  as  either  the  area  of  a 
diffracting  aperture  or  the  cross-sectional  area  of  a 
field  of  finite  extent.  The  far  field  approximation 
restricts  the  proximity  of  S  with  respect  to  S'  for  a 
given  area  A2  .  For  the  calculations  in  this  section, 
the  dimensions  of  S'  are  all  on  the  order  of  several 
microns,  the  wavelength  is  chosen  to  be  0.6  microns,  and  the 
distance  to  S  is  chosen  to  be  one  meter.  For  these 
parameters,  the  far  field  approximation  is  valid.  For  the 

__  A 

normally  incident  fields  considered  here  (k  *  k  az)  , 
most  of  the  energy  in  the  resultant  field  is  expected  to 
be  near  the  axis.  So,  we  also  make  the  axial  approximation 
that  6*0°  .  ’rhese  approximations  not  only  simplify 

the  calculations,  but  in  each  example  the  resultant  field 
is  then  proportional  to  either  a  Fourier  transform  or  a 
Fourier-Bessel  transform  of  the  initial  field.  For  the 
initial  fields  chosen  here,  the  transforms  are  solved  in 
closed  form  and  the  resultant  irradiance  patterns  are 
illustrated.  Finally,  the  Rayleigh-Sommerf eld  equation  is 


(59) 
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used  for  each  example  to  demonstrate  the  equivalence  discussed 
at  the  end  of  Section  III.  Let  us  then  begin  by  discussing 
the  linearly  polarized  field,  IT^ 

Example  1:  Linear  Polarization 

The  Solution.  Consider  the  field  defined  by 


on  S' 


everywhere  else 


C60) 


where  S'  is  the  rectangular  surface  defined  by  )x|  $  a 
and  |y|  <  b  .  Figure  3  shows  the  normalized  irradiance 
pattern  for  where  a  *  1.2  um  and  y  -  0  .  The 

magnetic  field  is  found  with  Maxwell's  equations.  Using 
the  definitions  of  Q  and  P  as  expressed  in  Eqs  (18) 
and  (19),  and  noting  that  specializing  for  S"  ,  z'  ■  0  , 

the  initial  field  is  represented  by 


% 


(61) 


P 


± 


(62) 


sin 


cos 


TTX  " 

7a- 


(63) 
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Figure  3.  Normalized  Irradiance  Pattern 

for  Initial  ^  with  a*1.2  ym, 

X-0.6  ym,  and  z»l  m. 
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where 


-  (k  ±  k  ) 
w  z 


(64) 


for 


k2 


kz2  *  <tt>2  *  C&>2 


and 


-•{ 


COS  cos 


Since  we  assumed  k 
and  the  ratio  of  F 


Vz 

to 


|x'|  5  a,  |y-|  !  b 
elsewhere 

«  2  ?r 2  r  1  1 

■  kl>> 

F+  can  be  expressed  as 


(65) 


FT 


x2  r  1  +  i i 


(66) 


For  a  i  A  and  b  >  A  ,  the  terms  containing  F_  can 
be  neglected  when  compared  to  terms  containing  F+ 
Neglecting  these  terms,  and  making  the  approximation 
discussed  earlier  that  e  ;  0°  ,  one  finds  that  only  the 

Q+  and  P_  integral  equations,  Eqs  (35)  and  (39),  are 
significant:  (See  Appendix  D) 
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Q  *  -P  -  -2  F.  k  /  E'  4>  dS' 

S' 


(67) 


Under  the  far  field  approximation  discussed  earlier,  the 
integral  can  be  written  as  two  Fourier  transforms  that 
have  closed  form  representations.  So  the  resultant  field 
is  expressed  as  Eq  (68). 


EjOO 


i  -  ikz 

17  e 


e 


ik 

Tz 


(x2+y2) 


TaB" 


cosfi^)  cos(i£) 

x  -  ax  (68) 

t(*)2  -  4)2k<*)2  -  t^)2] 

This  field  has  the  same  polarization  as  the  initial  field. 

When  the  arbitrary  assumption  is  made  that  the  components 
of  the  initial  field  can  be  treated  separately  with  the 
Rayleigh-Sommerfeld  approach,  the  same  resultant  field 
is  obtained  for  e  ■  0°  .  (See  Appendix  d) 

Utility  of  the  Fast  Fourier  Transform.  In  the 
derivation  of  Eq  (68),  the  Fourier  transform  in  Eq  (69) 
was  obtained  for  the  integral  of  x'  components  of  the 
initial  field. 
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«•  „  -  i 2  irf  X  * 

F(£x)  -  /  TCxOcos^-  e  x  dx'  (69) 

•  00 

where  T(x')  -  1  for  |x'|  5  a  and  is  zero  elsewhere,  and 
the  spatial  frequencies  are  f  ■  ~  .  The  integral  can 

a  A  2 

be  solved  in  closed  form  as 


FCfx) 


,  cost2'a  V 

2  a  2 


(70) 


An  alternate  approach  to  evaluating  the  integral  in  Eq  (69) 
is  by  numerical  computation.  Let  us  begin  by  reviewing 
some  essentials.  The  continuous  Fourier  transform  (CFT) 
pair  are  shown  in  Eqs  (71)  and  (72). 

«  i2irf  x 

f(x)  -  /  F(f  )  e  x  dfx  (71) 

•  - i2nf  x 

F(fx)  -  /  f(x)  e  x  dx  (72) 
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The  discrete  Fourier  transform  (DFT)  is  an  approximation 
to  the  CFT  that  allows  numerical  calculation  of  the  transforms. 
The  discrete  Fourier  transforms  pair  is  (Ref  17:41) 


N-l  lU  V 

f(k)  -  l  F  ( j )  e 

j-0 


.  N-l  -i27r  & 

F(j)  -  l  l  f(k)  e  N 

i,—  n 


,  k-0,1 ,. . . ,N-1 


,  j-0,1 , . . . ,N-1 


where  k  are  spatial  values,  j  are  spatial  frequencies 
and  N  is  the  number  of  sample  points  taken  of  the  function 
in  the  summation.  Direct  summation  of  the  DFT  requires  a 
time  proportional  to  N  ,  while  the  use  of  a  Fast  Fourier 
Transform  (FFT)  reduces  this  computation  to  a  time  propor¬ 
tional  to  N  log2  N  (Ref  18:70).  The  FFT  routine  used  for 
this  report  was  a  library  routine  in  the  International 
Mathematics  and  Statictics  Library  (IMSL). 

For  this  example,  the  resultant  field  was  calculated 
using  both  the  FFT  routine  and  the  analytic  expression 
which  permits  the  comparison  of  the  DFT  with  the  CFT.  In 
this  comparison,  an  important  tradeoff  is  made  between  the 
distance  over  which  the  initial  field  is  sampled  (the  data 
window)  and  the  number  of  sample  points  taken  in  this  window. 
Note  that  the  data  window  must  at  least  cover  the  non-zero 
portion  of  the  initial  field,  but  may  extend  beyond  this 
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region,  giving  a  finite  approximation  to  the  actual  geometry. 

In  the  discrete  case,  the  spatial  point  that  corresponds 

to  a  given’ spatial  frequency  f  is  given  by 

xi 

xi  «  (i)  (^-)  i*0 ,1 , . . .  ,N-1  (75) 

where  L  is  the  length  of  the  data  window  (Ref  19)..  So, 
for  a  given  N  ,  the  resolution  of  the  resultant  field, 
which  is  inversely  related  to  the  range  of  ,  increases 

as  L  increases.  But  the  comparison  for  this  example 
showed  that  (1)  for  fixed  N  ,  the  maximum  percent  difference 
between  the  CFT  and  DFT  values  increased  as  L  increased, 
and  (2)  for  fixed  L  ,  the  maximum  percent  difference 
decreased  with  increasing  N  .  One  concludes  that  the 
resolution  must  be  traded  against  an  accurate  representation 
when  using  the  DFT  to  compute  the  resultant  field. 

Figure  4  shows  the  normalized  irradiance  pattern  for 
the  resultant  field  using  the  FFT  routine.  Here,  N  *  128 
and  L  ■  12a  .  The  maximum  percent  difference  between 

this  result  and  the  normalized  irradiance  pattern  predicted 
by  Eq  (70),  shown  in  Figure  5,  is  1.08  percent.  The 
average  percent  difference  is  about  0.05  percent.  Thus, 
the  DFT  did  serve  as  an  accurate  approximation  to  the 
analytical  result  for  this  example. 
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Figure  4.  Normalized  Irradiance  Pattern  for  Resultant 
Using  the  Fast  Fourier  Transform. 


t.e 


Normalized  Irradi 


where  a  is  the  radius  of  a  circular  aperture  and 
kpa  *  3.831705970  (second  zero  of  Jj  (kpp)  ).  Figure  6 
shows  a  sketch  of  the  normalized  irradiance  for  an  aperture 
with  a  6  um  radius.  Using  E2  and  the  corresponding 
obtained  from  Maxwell's  equation,  Q'  and  P"  are  written 
on  the  surface  S'  .  These  can  be  written  in  component  form 
as  shown  in  Eqs  (77)  through  (79). 


<>;  - 

Ft e±iT'  w') 

(77) 

p;  - 

Fm  e*  lT'  J,(k  o') 

♦  1  P 

(78) 

ik  2 

v  - 

p;  -  -  — £-  J0(kpP') 

.(79) 

where  the  characteristic  scalar  is 
k 

F±  “  “  Ckz  1  k) 

u 

and 

k2  •  k2  ♦  k2  (81) 

z  p 


(80) 
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EXAHPLE  E2  (60URCE) 


Figure  6.  Normalized  Irradiance  Pattern  for  Initial  E. 

4 

with  a=6  pm,  A»0.6  pm,  and  z*lm. 


By  an  argument  similar  to  that  made  for  the  previous 
example,  terms  containing  F_  are  neglected  in  comparison 
with  terms  containing  F  +  when  a  is  on  the  order  of  a 
wavelength  or  greater.  By  using  the  six  scalar  integral 
equations,  Eqs  (35)  through  (40),  with  the  far  field 
approximation,  and  assuming  0  ;  0°  ;  the  resultant  field  is 


2na  -ikz 
\z  e 


jk£i  Jofk0a) 

CtE£T)!  -  « 

P 


(82) 


which  is  of  the  same  polarization  as  the  initial  field. 

The  Rayleigh-Sommerfeld  equation  yields  the  same  result  if 
one  assumes  the  transverse  rectangular  components  of  ^ 
can  be  treated  independently.  (See  Appendix  E)  Figure  7 
shows  the  resultant  field  for  a  -  6  ym. 


Example  3:  Circular  Polarization 

Owing  to  the  circular  nature  of  Q±  and  P  ,  this 
example  is  quite  straightforward.  Consider  the  right 
circularly  polarized  initial  field 

•ik  2  « ' 

r  e  a  for  p  <  a 

Ej  -  |  0  »  *  (83) 

0  for  p  >  a 
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EXWlfLE  E2  (OUTfUTl 


gure  7.  Normalized  Irradiance  Pattern  for  Initial 


EMHELE  E9  (SOURCE) 
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Figure  8.  Normalized  Irradiance  Pattern  for 
Initial  T%  with  a=6  ym,  A=0.6  yin. 


where  a  is  the  radius  of  the  aperture  and  kpa  *  2.40482S5577 
(first  zero  of  JQ(kpp)  )•  Figure  8  shows  the  normalized 
initial  irradiance  for  a  *  6  microns.  Then  on  S'  , 

Q'  ■  0  and  P'  *  0  exactly,  and 


q:  * 

T 

/T  i  F.  JQ(kpp') 

(84) 

p '  - 

i  F.  J0(IcpP'> 

(85) 

■ 

'  N 

O' 

i  k 

P'  ■  -  — —  (kpp')  (cos<>'  -  sin*') 

(86) 

where 

F*  " 

h  1  k«> 

u  z 

Again, 

k2  - 

k2  +  k2 
z  p 

(81) 

since  cylindrical  coordinates  are  used.  Recall  that  for 
+z  propagation,  P_  and  Q_  have  a  right  circular  nature 
as  reflected  here.  Again,  the  F+  terms  dominate  the  F_ 

terms  and,  using  the  six  integral  equations  (Eqns  (35)  through  j 

\ 

(40)}  with  the  far  field  approximations  and  letting  6  r  0°  ,  j 

the  resultant  field  is  i 
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ikp2 

2iri  a  J.  (k  a)  -v_  ”  Iz 

«  - i — e —  e  1KZ  e 

Azk 

P 


(87) 


which  is  a  right  circularly  polarized  field.  The  normalized 
irradiance  pattern  for  this  field  with  a  =  6  pm  is  shown 
in  Figure  9.  By  using  Eq  (54)  and  applying  the  Rayleigh- 
Sommerfeld  equation  to  each  transverse  rectangular  component, 
Eq  (87)  is  obtained.  (See  Appendix  F) 


Observations  on  Applications 

The  methodology  is  straightforward:  (1)  form  Q  and  P 
from  E  and  H  ,  (2)  use  the  component  equations  to  obtain 
the  resultant  Q  and  P"  ,  and  (3)  reform  E  and  R 
The  results  agree  with  scalar  diffraction  theory  in  the  far 
field  approximation  and  for  6  :  0°  ,  provided  the  scalar 

theory  is  applied  to  the  transverse  rectangular  components. 
Notice  that,  by  assuming  e  :  0°  ,  the  possibility  of  any 

off-axis  power  flow  is  eliminated  since  the  scalar  integrals 
that  could  contribute  to  such  power  flow  depend  directly 
on  9  .  The  possibility  of  polarization  shifts  for  the 

off-axis  region  is  similarly  negated.  An  example  of  such  a 

A 

polarization  shift  would  be  the  development  of  a  Sy  or 

A 

a  component  in  the  resultant  field  of  example  1.  Thus, 
the  complex  field  vector  approach  clearly  indicates  that 
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EXRHPLE  E3  (OUTPUT) 


Figure  9.  Normalized  Irradiance  Pattern  for  Resultant 


1 


scalar  theory  can  be  used  under  the  far  field  approximation 
and  for  e  ;  0° 

The  characteristic  scalar,  F  ,  not  only  simplifies 
the  application,  but  gives  insight  into  the  diffraction 
process.  If  the  aperture  is  smaller  than  two  wavelengths 
in  size,  the  F_  terms  could  not  be  neglected,  and  the 
resultant  field  would  be  changed  accordingly.  This  could 
give  rise  to  other  polarizations  and  change  the  amount  or 
direction  of  power  flow. 

Finally,  the  application  of  this  approach  is  straight¬ 
forward  since,  once  Q"  and  P'  are  formed,  the  resultant 
field  is  obtained  from  a  set  of  component  equations.  The 
application  is  particularly  straightforward  for  circular 
polarization  problems,  as  seen  in  example  3.  This  approach 
avoids  the  necessity  of  resolving  the  initial  E  into 
rectangular  components  as  required  when  using  scalar 
diffraction  theory.  However,  since  this  approach  agreed 
with  scalar  diffraction  theory  under  the  axial  and  far  field 
approximations,  and  since  the  scalar  theory  was  shorter  in 
application,  applications  in  this  region  are  more  quickly 
accomplished  with  scalar  theory  without  any  loss  of  accuracy. 
If  either  approximation  is  not  valid,  the  complex  field 
vector  approach,  though  more  lengthy  in  application,  is  more 
appropriate  than  the  scalar  theory.  So,  this  approacl;  ’is 
not  only  rigorous  in  its  development  as  shown  in  Section  III, 
but  it  is  applied  in  a  straightforward  manner  as  demonstrated 
in  this  section. 


47 


Conclusions  and  Recommendations 


In  this  section,  three  important  conclusions  and  four 
major  recommendations  are  discussed.  As  a  brief  summary, 
recall  that  the  use  of  the  complex  field  vectors  led  to 
uncoupled  field  equations.  Then  these  equations  were  solved 
by  rewriting  the  results  obtained  by  Doughty  in  terms  of 
Q  and  P  .  The  resultant  field  was  expressed  at  this 
point  by  a  pair  of  uncoupled  vector  integrals  over  the  open 
surface  S'  .  When  S'  was  planar,  the  resultant  field 
was  written  in  terms  of  six  component  integrals  over 
S'  .  Section  IV  considered  three  examples  that  showed  how 
easily  these  integrals  are  applied.  In  the  special 

a 

case  where  (1)  k  Z  k  a  initially,  and  (2)  0  <<  1  rad  , 

the  tangential  components  of  the  resultant  field  were 
expressed  in  a  form  identical  to  the  Rayleigh- Sommerf eld 
equation.  With  this  summary  in  mind,  the  following  conclu¬ 
sions  are  presented. 

Conclusions 

First,  the  objectives  set  forth  in  Section  I  have  been 
met  in  this  report.  The  complex  field  vector  approach  leads 
to  two  separated  vector  integral  equations,  Eqs  (24)  and  (2S) 
that  are  evaluated  over  an  open  surface  S'  .  The  only 
restrictions  on  these  equations  are  the  assumptions  that 
Ohm's  Law  holds  in  the  medium  and  that  the  medium  is  homogen¬ 
eous,  isotropic,  linear  and  free  of  charge,  and  the 


requirement  that  the  initial  E  and  H  must  satisfy 
Maxwell's  equations  on  S'  .  When  S'  is  specialized  to 
be  a  plane  surface,  the  six  component  integrals  in  Eqs  (35) 
through  (40)  were  derived  by  the  further  restriction  that 
k  >>  1/R  .  Thus,  the  development  of  resultant  field  is 

rigorous,  as  is  required  for  an  acceptable  vector  diffraction 
theory.  This  approach  has  the  further  advantage  that  the 
application  is  mathematically  simple,  since  only  separated 
component  integrals  are  involved.  Indeed,  for  the  special 
case  where  k  :  k  az  initially  and  e  <<  1  rad  ,  the  scalar 
equations  express  the  transverse  components  of  the  resultant 
E  in  a  form  identical  to  the  Rayleigh-Sommerfeld  equation. 

So  the  complex  field  vector  approach  can  be  as  simple  and 
straightforward  as  scalar  diffraction  theory. 

Second,  this  derivation  of  the  Rayleigh-Sommerfeld 
equation  explains  why  this  equation  accurately  describes  the 
near-axial  irradiance  pattern.  It  further  implies  that  the 
Rayleigh-Sommerfeld  equation  is  not  accurate  in  the  region 
farther  off  the  axis.  In  this  region,  the  set  of  scalar  integral 
equations  should  be  used.  Thus,  the  complex  field  vector 
approach  is  clearly  more  general  than  scalar  diffraction  theory. 

Third,  the  applications  discussed  in  Section  IV  provide 
insight  into  the  diffraction-  process.  The  characteristic 
scalar  F  that  occurred  in  the  three  examples  describes  the 
effect  of  the  aperture  size  on  when  terms  in  F_  could  be 
neglected  by  comparison  to  terms  in  F+  .  Also,  the 
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circularly  polarized  nature  of  Q+  and  P+  that  was 
discussed  in  Section  III  is  strongly  evident  in  the  third 
example  which  dealt  with  a  circularly  polarized  initial 
field.  Further,  it  is  important  to  note  that  the  terms 
that  could  have  caused  polarization  shifts  and  off-axis 
power  flow  were  eliminated  only  by  the  angular  restrictions 
and  not  inherently  by  the  component  integrals.  Finally,  the 
restriction  that  e  ;  0°  removed  the  terms  containing  & 
(Recall  that  6  ■  tan”1  .  The  ratio  of  coordinates 

relates  field  values  at  a  point  on  the  initial  field  to  the 
value  of  the  field  at  an  observation  point,  which  gives  a 
certain  orientation  between  initial  and  resultant  fields.) 

So  it  is  not  surprising  that  the  results  agree  exactly 
with  the  predictions  of  scalar  diffraction  theory. 

Based  on  these  three  main  conclusions,  the  following 
recommendations  are  made. 

Recommendations 

First,  the  calculation  of  off-axis  power  flow  and 
polarization  shifts  needs  to  be  developed  more  fully. 

Using  the  Fresnel  approximation  in  place  of  the  far  field 
approximation  casts  the  component  integrals  in  the  form  of 
two  dimensional  convolutions.  These  equations  may  be 
amenable  to  fast  Fourier  transform  techniques  if  the  initial 
field  is  chosen  properly.  The  resulting  components  of  $ 
and  7  would  then  be  used  to  form  7  and  IT  as  well  as 
the  Poynting  vector  7-7x7  .  From  these  three 


50 


quantities,  the  time  average  polarizations  and  power  flow 
can  be  studied. 

Second,  the  sensitivity  of  this  approach  to  the 
requirement  that  the  initial  E  and  H  satisfy  Maxwell's 
equations  on  S'  should  be  assessed.  For  example,  if 
only  the  initial  transverse  components  of  E  and  the 
direction  of  propagation  were  known,  some  information  about 
the  resultant  field  may  still  be  obtained  by  partially 
forming  Q  and  P  .  A  weakening  of  the  above  requirement 
would  make  the  results  more  useful  in  practical  applications. 

Third,  an  example  should  be  worked  where  the  field  is 
discontinuous  at  the  edge  of  the  surface  S'  ,  so  that  an 
edge  current  is  introduced  in  order  to  satisfy  Maxwell's 
equations  on  S'  .  The  example  would  illustrate  how  to 
deal  with  these  fields.  The  resultant  field  should  be 
compared  against  the  field  predicted  by  the  Rayleigh- 
Sommerfeld  equation  to  assess  the  effect  of  edge  currents 
on  diffraction. 

Fourth,  the  term  6  discussed  in  the  third  conclusion 
should  be  investigated.  If  the  initial  field  is  azimuthally 
symmetric,  «  would  seem  less  significant  since  the  spatial 
phase  shift  should  also  be  azimuthally  symmetric.  Also,  the 
ratio  of  coordinates  „  is  not  unique  for  each  set  of 
points  described  by  r-r '  .  This  mapping  should  be  investi¬ 

gated  to  develop  more  clearly  the  nature  of  6 

Throughout  this  report  it  has  been  emphasized  that  the 
diffraction  problem  has  been  solved  by  the  complex  field 
vector  approach  in  a  manner  that  is  both  rigorous  and 


mathematically  simple.  Other  solutions  to  the  diffraction 
problem  exist,  some  more  useful  than  others.  Unfortunately, 
the  most  useful  (scalar  diffraction  theory)  is  the  least 
rigorous  and  the  more  rigorous  approaches  are  sufficiently 
complex  as  to  be  less  useful.  The  complex  field  vector 
approach  spans  the  "middle  ground"  of  being  derived  directly 
from  Maxwell’s  equations  and  being  useful  in  obtaining  both 
solutions  and  insight  into  diffraction  problems. 
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Appendix  A 


Components  of  the  Vector  Integral  Equation  for  (} 

The  purpose  of  this  appendix  is  to  obtain  Eq  (26)  from 
Eq  (24).  Let  us  begin  by  recalling  Eq  (24). 

Q(r)  *  /  {[nxQ'(r)] *vxG  ♦  r-[nxQ' (r ')  ]  •  vxV'xG-}  dS'  (A.l) 

S'  K 

Letting  nxQ'(r')  ■  J(r')  ,  this  equation  can  be  written  as 


Q  -  /  {vx[J*C  +  i-  J'V'xC]>  dS' 

S'  K 

Recalling  that  G  *  T*  ,  one  readily  finds  that  v'xC 
is  the  completely  asymmetric  dyad. 


(A.  2) 


*\  A 

where  ♦  =  rr  ,  etc.  Then  by  direct  manipulation, 

)  •  0  Z 

one  has 


S5 


Vx(J-7'xf)  -  ax  [Jy7,xy  *  V,xz  -  Jx(*>yy  *  ♦  _„)] 


+  a  [J  $  +  J  #  -  J  (♦  ♦  $  )1 

y  L  z  ,yz  x  ,yx  yv  ,xx  ,zz'J 


♦  a  [J  *  +  J  *  -  J  (♦  ♦  *  )]  (A. 4) 

z  L  x  ,zx  y  ,zy  z  ,xx  ,yy  J  v 


Observing  that  J*?T  *  J*  ,  the  other  term  in  the  integrand 
is  written  in  component  form  as 


VXJ«  -  ax  (J2.>y  -  Jy.<z) 


*  ay  <Jx\,  '  V,x> 


+  a  (J  *  "J#) 

z  y  ,x  xv,y ^ 


(A.  5) 


Noting  that  *  satisfies  the  scalar  Helmholtz  equation, 

*  A 

Eq  (14),  and  choosing  n  ■  +a  so  that  S'  is  planar  and 


J  "  -<Vax  +  Qxay 


(A. 6) 


one  obtains  for  the  resultant  field 
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I  <-y  j 

S'  *  * 


/  <-<C« 


S' 


y  .2 


I  (Q't  x 

g  /  A  >  A 


-  ^Cq^C^2*+*(XX)  -  Q'»>xyl>  ds 

-  jr[-Q'(k2t*«>yy)  *  «' 

*  V.y  *  JF[V,2y-y.2*])  dS' 


(A. 7) 


A 

Taking  the  scalar  product  of  ax  with  Eq  (A. 7),  Eq  (A. 8) 
is  obtained: 


1' 


<-v,* 


-  ^-CQyCk2**^ 


xx 


)  - 


^x®,xy^ 


dS' 


(A. 8) 


This  is  Eq  (26)  as  desired.  Equations  (27)  and  (28)  are 

A 

obtained  by  the  scalar  product  of  Eq  (A.  7)  with  a^,  and 

A 

az  ,  respectively. 
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Appendix  B 

Derivatives  of  the  Green's  Function 

This  appendix  obtains  the  first  and  second  derivatives 
of  the  scalar  portion  of  the  dyadic  Green's  function, 

*  *  exp  [-ikR]/4irR  .  The  first  derivative  with  respect 
to  x^  is  obtained  directly  as 

(x .  ■  x') 

®,x.  "  f£.  *  -  (ik  +  1/R)*  (B.l) 


where  x ^  is  x  for  i  *  1  ,  y  for  i  *=  2  ,  and 
z  for  i  -  3  .  The  second  derivative  is  more  complicated. 
By  the  rules  for  differentiating  a  product,  one  has 


*,xixj  *  wr  ^ 


(ik 


1 


3(xi-xf) 

a3T 

J 


(x^xr)  (Xj-xp 
R2 


i  2  (x^-xpcx.-xo 

-  (ik  ♦  Jo2*  -1—1 - 1— L- 

K  R2 


(x.-x() (x.-x() 

♦  i  .  .  -  J.  l —  $ 

R3 


(B.  2) 


3(x.-xp 

Noting  that  — - -  -  6..  ,  and  by  direct  algebraic 

0X4  1J 


manipulation. 


3 
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6.  .  ,  (x •  -X?)  (X.-X.'j  ,  , 

*»xixj  "  +  R-)  ♦  —  -  -2— — 2-[ik+£(ik*l)-k2+i2]}  (B.  3) 


Under  the  approximation  k  1/R  ,  Eq  (B.l)  becomes 


-  ik  4 


(x,-xr) 


(B.4) 


and  Eq  (B.3)  becomes 


■xixj 


6..  (x.-xp(x.-xj) 

-  ik  *  -ii-  -  k2  — i — — 2— 

R  n  2 


(B.S) 


These  last  two  equations  are  Eqs  (29)  and  (30). 


Appendix  C 

Expression  of  Complex  Component  Q+  Using  Impedances 
We  have  the  complex  component 


Q  ■  Q  +  i  Q 

X+  Xx  Xy 

where 

Q  -  wH  ♦  i/ul  E 
x  x  * 

and 

Qy  ■  v*Hy  +  i/7F  Ey 


Substituting  Eqs  (C.2)  and  (C.3)  into  Eq  (C.l),  and 
factoring  out  the  component  of  E  , 

q4  -  ©Ej  c-1  ♦  i m  ex  a  *^/f 

Using  a  complex  form  for  the  intrinsic  impedance, 


CC.l) 


(C.2) 


(C.3) 


(C.4) 


(C.  5) 


and  defining  the  wave  impedances  following  Harrington, 

I 


(Ref  16:86) 
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where  the  upper  sign  is  for  +z  propagation  and  the  lower 
sign  for  -z  propagation,  one  obtains  the  final  expression 
for  Q+  that  is  Eq  (49): 


Q+  -  -/pT  Ey  (1  )  ♦  i/jiT  (1±£  )  (C.8) 

yx*  xy* 


Appendix  D 

Development  of  Example  1 


This  appendix  gives  a  more  detailed  development  of 
Example  1,  where  the  field  was  defined  by  Eq  (60).  On 
S'  ,  one  has 


H  "  cos  ft  cos  ft  ax 

R1  *  ay(kzcosft  cosft) 

Using  the  definitions  of  Q  and 
Eq  (34),  the  components  of  Q 


(D.l) 

+  iaz(ft  sinft  cosft)  (D*2) 

P  ,  Eq  (18)  and  (19),  and 
and  P'  are 


- 

1  F±  Ex 

■  (D.3) 

p ;  - 

-i  F  E' 

(D.4) 

♦ 

- 

n  /  1  IT  •  TTY  irX 

Pz  TuF  sin  ft  cos  7a 

(D.S) 

where 

F±  ■  h  <k  1  (»•« 

and 
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cos 


(D.7) 


E; 


1TX  ' 

la 


cos 


r 


As  discussed  in  Section  IV,  the  F_  terms  can  be  neglected 
when  compared  to  the  F+  terms  for  a  >  X  .  Substituting 
the  value  in  Eq  (D.3)  into  the  scalar  integral  equations, 

Eqs  (35)  through  (37),  and  making  the  paraxial  approximation 
of  Eq  (56),  one  finds 

(D.8) 

(D.9) 

(D. 10) 

A  similar  substitution  is  made  using  Eq  (D.4)  for  the  P 
equations.  One  should  note  that  all  the  components  of  Q 
and  P  are  present  in  the  resultant  field,  so  all  the 
components  of  E  and  H  are  possible.  This  allows 
polarization  shifts  and  off  axis  power  flow  to  occur.  Also, 
note  that  all  the  terms  containing  6  are  proportional 
to  e  ,  so  as  e  approaches  zero,  these  terms  vanish  and 
the  sense  of  orientation  represented  by  5  is  lost.  Let 
6  Z  0°  Then,  for  Q 


Q+  :  -k  /  F.  (2-02)  E-  *  dS' 

S'  * 

Q  ;  -k  /  F  y2eZl6  E'  •  dS' 

S'  x 

<JZ  =  -  -r  LF-  eis  6  E; » 

O 
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<*♦ 


Q. 


F  k  /  E'  «  dS' 


and  for  F  , 


p-  :  -Q+ 


p+  :  pz  :  o 


The  integral  in  Eq  (D.ll)  is 


b  a  ,  , 

i  -  j^e'  «  dS"  -  [bdy'{  dx'  cosll  cos  S  * 


Under  the  far  field  approximation,  Eq  (59),  and  since 
0  ;  0°  ,  the  *  term  is  approximated  by 


e"ikR  .  l  -ikz  ’  ^x2+y2)  ~(xx'*yy-) 

T 7JT  *  777  e  e  e 


So  the  integral  becomes 


ik 


I  .  «'ikJ  zz(x2*y2)  Ix  I, 
1  “77 T  e  x  > 


where 


(D.ll) 

(D. 12) 


(D. 13) 

(D. 14) 


(D. IS) 


(D, 16) 


(D. 17) 
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a  • kxx ' 

Iv  ■  /  cos  rj  e1  z  dx' 


and  similarly  for  I 


Consider  the  integral  I  .  Let  x'  be  -x'  and 
multiply  Eq  (D.18)  by  an  aperture  function  T(x')  where 


T(x')  - 


1  |xi  f  a 


0  lx*)  >  a 


Then 


-  /  [T(x')  cos^f]  e  z  dx‘ 


This  is  the  Fourier  transform  of  [T(x')  cos^j]  , 


Ix  -  F{T (x')cos  } 


which  has  a  closed  form  solution 


,  o(fe&, 

11  cc#2  -  #)h 


(D.18) 


(D. 19) 


(D. 20) 


(D-21) 


A  similar  analysis  for 


yields 


IT 
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cos 


(i£> 


ct*)2  - 


So  the  resultant  field  is  expressed  by 


ik, 


iv  -ite  TT'1'2*)'2’  .  . 

■TF*  e  J-  1 


x  y  x 


where  Eq  (41)  was  used.  Noting  that  F+  ;  —  ,  one 


finds 


w„_  i  -ikz(x2+y2)  tt2 

E  -  ”e  *55 


cos  fecosfM) 


\z 


cc^)2  -  c^)2][(^)2  -  4)2] 


This  is  Eq  (68). 

If  one  assumes  the  transverse  components  of  the 
initial  field  can  be  treated  independently  with  the  Rayleigh- 
Sommerfeld  equation,  Eq  (17),  then  the  resultant  field  is, 
for  0  Z  0°  and  using  the  far  field  approximation. 


i  -go-V2) 

Tz 


*  =  --e  «'  "  /  cosgcosfe 


.  i|(xx'+yy-) 


dS 


(D.23) 


(D.  24) 


(D.25) 


(D.  26) 


This  integral  is  the  same  one  evaluated  in  Eq  (D.1S),  so 
the  resultant  field  obtained  from  scalar  diffraction  theory 
is  the  same  as  that  obtained  in  Eq  (D. 25). 


Appendix  E 

Development  of  Example  2 

This  appendix  gives  more  details  about  Example  2, 
where  the  field  on  S'  was 


*2  “ 

-k  J.(k  p')  a 

pip 

k  k  ik  2 

A 

(E.l) 

«2  * 

-£-5-  J.  (k  p')  a  -  — J n(k  P') 
up  1  p  up  0  p 

az 

(E.2) 

One  forms  the  components  of  Q'  and  P'  by 

definitions : 

using  the 

F_  e*1*'  J,  (k  p') 

±  l  p 

(E.3) 

- 

F  e*1*'  J,  (k  p') 

;  1  P 

ik  2 

(E.4) 

v  - 

where 

Pz  *  ' 

k 

(E.  5) 

F*  " 

-f  (K  ±  k) 

u  z 

(E.6) 

As  argued  in  Section  IV,  the  F_  terms  are  neglected  in 
comparison  to  F+  terms  so  long  as  a  i  \  .  The  terms 

in  Eqs  (E.3)  and  (E.4)  are  substituted  into  Eqs  (35)  through 
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(40).  The  discussion  in  Appendix  D  applies  here  for  the 
paraxial  approximation.  Let  us  proceed  immediately  by 
letting  0  :  0°  .  The  far  field  approximation  for  *  in 


cylindrical  coordinates  is 


« 


k  p' 

+i-|-  cos(t'-t) 


(E.7) 


So  the  resultant  field  is  described  by  the  components 


Q+ 

2ikF  + 

4irZ 

-  ikz 
e 

ik  2 

P 

e  2z  I+ 

(E.8) 

2ikF+ 

-  ikz 

ik  2 

p 

(E.y) 

P_ 

4irZ 

e 

e  I_ 

Q_ 

Z  Q  Z 

Hz 

p+  : 

O 

1  • 

M 

P* 

(E. 10) 
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where  the  identities 


J.„(x)  -  (-l)n  Jn(x) 

Jn(-x)  -  C-l)n  Jn(x) 

were  also  useful.  With  the  use  of  an  aperture  function 
like  Eq  (D.19),  the  integral  in  p'  is  a  Fourier-Bessel 
transform  (or  a  Hankel  transform  of  order  one)  that  may  be 
solved  in  closed  form.  So  the  resultant  field  is 


ik  2 


-a 


2ira 

\z 


-  ikz 
e  e 


kap  ■ 


‘  Tz*  JlC"T")J0(kpa) 


p 


-  1 


where  the  relation  F+  Z  was  used.  This  is  Eq  (82). 

In  order  to  apply  scalar  diffraction  theory,  one  must 
first  write  the  initial  field  in  rectangular  coordinates, 
then  the  initial  field  on  S'  is 

Ez(r')  -  ax[kpJ1  (kpP')sin^']  -  ay [kpJx (kpp ')cos* '] 

Then,  assume  each  component  can  be  treated  independently 
with  the  Rayleigh-Sommerfeld  equation,  Eq  (17),  and  use 
the  exponential  forms  of  the  circular  functions. 


(E.14) 

(E.  IS) 


(E. 16) 


(E* 17) 
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sin  $ 


cos  $ 


(E. 18) 


-i* ' 


e^'-e"^ 


Consider  the  ax  component  first.  Under  the  far  field 

A 

approximation  and  letting  e  :  0*  ,  the  resultant  ax 

component  is 


E 


x 


i  3 


(E. 19) 


where  the  integrals  I  are  those  evaluated  above  in 


Eq  (E.ll).  The  resultant  Ex  is  then 

ik  2 

C  _  2ira  -ikz  *'ir  T  a^cin.  “1 ' 


j,  fifi) 


A  Z 


(ft)  -  1 

p 


(E. 20) 


From  a  similar  analysis,  E^  can  be  found  and  one  recombines 
the  components  using  the  relation 


E.  -  -E  sin$  +  E„cos«ji  (E . 21) 

<p  x  y 

and  the  resultant  field  in  Eq  (E.16)  is  again  obtained.  One 
should  note  the  necessary  step  of  casting  the  initial  field 
in  terms  of  rectangular  components  before  the  Rayleigh- 
Sommerfeld  equation  can  be  correctly  applied.  This  is  due 

A  A 

to  the  fact  that  a  and  a  can  be  factored  out  of  integrals. 

X  / 
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Appendix  F 

Development  of  Example  3 

In  this  appendix.  Example  3  is  developed  in  more 
detail.  The  initial  field  was  described  by  Eq  (83).  On 
S  ,  we  have  from  Maxwell's  equations  and  Eq  (54), 


E'  -  !  Jfl(k  o')a  *  -  J0(k  p')av 

sz  0  p  x  n  0  p  y 


(F.l) 


-  JL  Ik  j  (k  p')av  +  -L  !k  j  (k  p ")a 

wu  yy  O  p  x  wy  O  p  y 


\  k  * 

- 2-  J.  (k  p')(cos*'  -  sin*')a 

uUyjlp  T  '  z 


(F.2) 


The  components  of  $  and  7  follow  directly: 


q;  - 

0 

(F .  3) 

q:  * 

n  i  F.  J0(kpp') 

(F.4) 

p;  - 

0 

(F.  5) 

p  *  « 

-/I  i  F,  J0(kpp') 

CF.6) 

Q-  - 

ik 

p;  ■  -  — £j. (k  p ')  (cos*'  -  sin*') 

2  /T  u>  p 

(F.  7) 
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where  F±  ■  i-  (k  ±  kz)  .  As  argued  in  Section  IV,  the 
F_  terms  are  neglected  in  comparison  to  the  F+  terms. 
Then,  substituting  into  the  scalar  integral  equations, 
Eqs  (35)  through  (40),  with  the  far  field  approximation 
and  letting  e  :  0°  ,  one  finds  Q  =  0  ,  P  «  P  »  0 

and 

ik  2 

P.  :  -  A)  e-ikI  e"7*"  I 


where 

kpp  ' 

a  Zv  i — -—cos  ($'■$) 

I  -  /  J0(k  P')P'dp'  /  e  2  d*' 

o  p  o 

Using  the  identity  shown  in  Eq  (E.12),  the  integral  becomes 

a  k  p  ' 

I  -  -2*  /  p'J0(k  P')J0(-%-)  dp' 

o 

With  the  use  of  an  aperture  function  like  that  described 

in  Eq  (D.19),  this  integral  is  a  Fourier- Bessel  transform 

(or  a  Hankel  transform  of  order  zero)  of  Jrt(k  p')  •  The 

u  p 

integral  has  a  closed  form  solution,  so 


4_/2jra 

IzTv 


J1  (kpa) 


■ikz 


k  2 

.  i  p 

1_ZT 


k  a 

v-H 


r  y*  i2 

in5 

p 


(F.8) 


(F.9) 


.10) 


(F.ll) 


72 


From  Eqs  (41)  and  (42),  one  finds 


m 

so  the  resultant  field  is 


(F. 12) 


J3  * 


2m  a  Jj(lc  a)  ik2  7§- 

- ITT -  c  e 


k  a 

Vt"> 


/•  kp  \ 2 

■tnr0 


-  i 


(F.  13) 


This  is  Eq  (87). 

To  use  scalar  diffraction  theory,  one  must  first 
obtain  the  rectangular  components  of  E  as  in  Eq  (F.l), 
and  then  assume  each  component  can  be  treated  independently. 
The  equation  for  Ex  is  (for  9  Z  0°  ) 


-at  |.  W'> 


■ikr 

U — 


dS- 


(F. 14) 


With  the  far  field  approximation,  one  has  the  same  integral 
that  was  solved  earlier  in  this  appendix.  The  resultant 


/r  irk  ai 

p 

xz 


0 


Jl(kpa) 


ik  2 

_ p _ 

ikz  ^  2z 

e  [(l£)* 


k  a 

c4- 


) 


(F. IS) 
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A  similar  development  for  E^  leads  to 


nk  a 


Az 


£ —  J,  (k  a)  e'ikz  e 
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Then,  using  the  relation  E  *  —  (E  -  iE  )  ,  Eq  (F.13) 

/I  ' 

is  again  obtained  for  the  resultant  field.  Thus,  scalar 
diffraction  theory  and  the  complex  field  vector  approach 
yield  the  same  result  in  this  region. 
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